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Abstract 

We calculate the unpolarized QED Compton scattering cross section in a manifestly covariant 
way. Our approach allows a direct implementation of the specific kinematical cuts imposed in the 
experiments, e. g. HERA-HI. We compare the 'exact' cross section in terms of the structure func- 
tions Fi^{xb, Q 2 ), assuming the Callan-Gross relation, with the one obtained using the equivalent 
photon approximation (EPA) as well as with the experimental results. We find that the agreement 
with the EPA is better in x 7 bins, where x 7 is the fraction of the longitudinal momentum of the 
proton carried by the virtual photon, compared to the bins in the leptonic variable xi. 



*Electronic address: asmita@physik.uni-dortmund.de 
^Electronic address: pisano@physik.uni-dortmund.de 



1 



I. INTRODUCTION 



The QED Compton scattering in high energy electron-proton collisions e(l)p(P) — > 
e{l') r y{k')X, with a real photon j(k') emitted at the lepton vertex (Fig.l), is one of the 
most important processes for an understanding of the photon content of the proton in the 
so-called 'equivalent photon' approximation (EPA), first introduced for a charged particle 



Weizsacker and Williams 
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l| and later extended to the nucleon and investigated widely 



. QED Compton scattering has been recently analyzed by the HI collab- 
oration at HERA j^. In addition to the information concerning the photon content of the 
proton within the framework of EPA, it can also shed some light on the structure functions 
Fj^(x B , Q 2 ) of the proton P, Q] in the low-Q 2 region, which are presently poorly known 
In |9| these alternative descriptions were confronted with data and it was found that 



0. 



the description in terms of F 12 , i.e. for X ^ p, is superior to the one in terms of the 
inelastic photon distribution 7inei(^_B ) Q 2 )- Henceforth we shall refer to the description in 
terms of as 'exact' to distinguish it from th e ap proximations involved in the EPA. It 
should be noted, however, that the analysis in l| [l(| utilized the Callan-Gross relation l^J 
Fl{xb, Q 2 ) = F 2 (xb, Q 2 ) — 2xbF\(xb, Q 2 ) = 0. This relation is contaminated by higher or- 
der (NLO) QCD corrections as well as by higher twist contributions relevant in the low-Q 2 
region which may invalidate the assumptions underlying the 'exact' analysis. Furthermore, 
the analysis in jjj, Q| was carried out within the framework of the helicity amplitude formal- 
ism [10]. The implementation of experimental cuts within this formalism is nontrivial and 
affords therefore an iterative numerical approximation procedure jloL Q] whose first step 
corresponds to —k 2 = Q 2 = 0, where k is the momentum of the virtual photon. 

It is this second issue which we intend to study here. We shall replace the noncovari- 
ant helicity amplitude analysis of [10] by a standard covariant tensor analysis whose main 
advantage, besides compactness and transparency, is the possibility to implement the exper- 
imental cuts directly and thus avoid the necessity of employing an iterative approximation 
of limited accuracy. The first issue concerning the Fl contributions affords some estimates 
of this poorly known structure function and we refrain from its study here. 

The plan of the paper is as follows. In section II, we calculate the exact cross section for 
the elastic scattering. In section III, we calculate the cross section for the inelastic channel. 
Our numerical results are discussed in section IV. The summary is given in section V. All 
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the useful and necessary formulae and the kinematics are given in the appendices A, B, C 
and D. 



II. ELASTIC QED COMPTON SCATTERING 

We consider elastic QED Compton scattering: 

e(l)+p(P)^e(l') +1 (k>)+p(P'), (2.1) 

where the 4-momenta of the particles are given in the brackets. We introduce the invariants 

S = (P + l) 2 } s = (l + k)\ t = k 2 . (2.2) 

Here k = P — P' is the 4-momentum of the virtual photon. The photon in the final state 
is real, k' 2 = 0. We neglect the electron mass everywhere except when it is necessary to 
avoid divergences in the formulae and take the proton to be massive, P 2 = P' 2 = m 2 . The 
relevant Feynman diagrams for this process are shown in Fig. 1, with X being a proton and 
Px = P' ■ The squared matrix element can be written as 

I M el | 2 = ~flS"(P, P')T, U {1, k; I', k'), (2.3) 

where 

H^(P,P') = \Y,( P ' I ^(0) I P)*(P' I AO) I P) (2-4) 

spins 

is the hadronic tensor, J M being the electromagnetic current. 
If we use the notation 

N N d 3 v- 

dPs N (p; Pl ,..., PN ) = (27r) 4 tf(p-x; ft )n (2^%p ( 2 - 5 ) 

for the Lorentz invariant A^-particle phase-space element, the total cross section will be 

^(S) = 2{s \ m2) J dPS 2+1 (l + P; /', k', P')\M el \ 2 . (2.6) 

Eq. ()2.6j) can be rewritten following the technique of Q], which we slightly modify to imple- 
ment the experimental cuts and constraints; in particular all the integrations are performed 
numerically. Rearranging the (2 + l)-particle phase space into a sequence of two 2-particle 
ones, Eq. ()2.6|) becomes: 

1 r ds 1 
<ra{S) = 2{s _ m2) J ^ dPS a (l + P; r + k', P')^{P, P')X,„(l, k) . (2.7) 
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X^v contains all the informations about the leptonic part of the process and is defined as 



X^l, k)=J dPS 2 (l + k- r, k')T^ u (l, k- 1', k'), 



being the leptonic tensor [1J, [15J] : 



Tpi/il, k]l,k) 



4e 4 f 1 



su I 2 



<v (s 2 + u 2 + 2«) + 2s z M z„ + 2m z;x 
+ (t + t)(ix + z„zy - (5 - 1) (zx + z^y 

+ (« - 1) («x + tk') 



(2.8) 



(2.9) 



where we have defined t = (/ — Z') 2 and u = (I — A;') 2 . It can be shown that 

dPS 2 (l + k;l' } k')= dt 2 f V * (2.10) 
16 7r(s — t) 

with yj* denoting the azimuthal angle of the outgoing e — 7 system in the e — 7 cm. frame. 
For unpolarized scattering, X^ v is symmetric in the indices //, z/ and can be expressed in 
terms of the two Lorentz scalars X% and X 2 : 

X»v(l,k) = I|[3X 1 (a,t)+X 2 (S l t)](^i-Ar) (j^Z - ^ 



+ [Xi(s,t) +X 2 (s,t)](%,-^^) 



(2.11) 



with 



4t 



(S-t) 



7 Frx^ u (i,k), 



(2.12) 



X 2 (s,t) =g^X^ u (l,k). 

Using the leptonic tensor (|2.9j) and also the relations 

l-k = -(s-t), I ■ P = -{S - m 2 ), k-P = -t, 
2 V 1 2 y ' 2 

we obtain 



(2.13) 



(2.14) 



t l^l Tui/ A tt 

- e 



4vr 2 (s-t) 3 2n 2 (s-tf 



x 1 (§,t,i), 



(2.15) 



g^T^ 4 (t 2 - 2ts + 2s 2 + 2st + i 2 ) 



16vr 2 (s - t) 



4tt 2 s(s - t)(t - s - t) 



x 2 {s,t,i), 



(2.16) 
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where e 2 = 4ira. The invariants Xi(s,t,i), with i — 1,2, are related to X;(s,t) by 

Xi(s,t)=2irl dtXi(s,t,t). (2.17) 



ft 

' ^min 

The integration limits of i are: 



tmax = 0, t min = -s + t + — - m 2 , (2.18) 

where m e is the mass of the electron. We point out that the kinematical cuts employed by us 
prevent the electron propagators to become too small and thus the divergences are avoided, 
so we can safely neglect the electron mass in the numerial calculation. The hadronic tensor 
in the case of elastic scattering can be expressed in terms of the common proton form factors 
as 

H^(P, P') = e 2 [H 1 (t)(2P - ky(2P - kf + H 2 (t)(tg^ - k"k v )\, (2.19) 

with 



Hl{t ^m T m^3M, H 2(t) ^ M(t) . ( 2 ,o, 



The electric and magnetic form factors are empirically parametrized as dipoles: 



Using 



Gfi(f)= [l- f /(0.71GeV')r <*,<*)= 2.79 G,(t). (2.21) 



dPS 2 (l +P;V + k', P) = - f (2.22) 



finally we get 



a 



°cl(S) = — ^ ^ / ds — I dt dip" 



S -m 2 fS -m 2 

2— — ; 1 

s — t V s — t 



2m 



2 



x [3Xi(s, t, t) + X 2 (s, t, t}} + — [X^s, t, t) + X 2 (s, t, t)} + X 1 (s, t, t) 



Hi(t) 

+ X 2 (s,t,t)H 2 (t)\, (2.23) 



where s m i n denotes the minimum of s and t niint max are given by 

1 



t ■ — 9m 2 — 

l mm, max ^ui 



(S + m 2 )(S-s + m 2 ) ±(S- m 2 )^J{S - s + m 2 ) 2 - ASm 2 



[2.2A) 



It is to be noted that in Eq. (J2.23JI we have shown the integration over ip* explicitly, because 
of the cuts that we shall impose on the integration variables for the numerical calculation 
of the cross section. The cuts are discussed in section IV. The EPA consists of considering 
the exchanged photon as real, so it is particularly good for the elastic process in which the 
virtuality of the photon \t\ is constrained to be small (< 1 GeV 2 ) by the form factors. It is 
possible to get the approximated cross section c^j^^ from the exact one in a straightforward 
way, following again |2| . If the invariant mass of the system e — 7 is large compared to the 
proton mass, s min ^> m 2 , one can neglect \t\ versus s, m 2 versus S, then 

Xl(s, t, i) « Xi(s, 0, t) = 0, (2.25) 

and 

X 2 (s,t,t) *X 2 (s,0,t) = -H^lfi, (2.26) 

7r at 

where we have introduced the differential cross-section for the real photoproduction process 
e7 — > ej: 

dais, i) 2na 2 ( u s\ , . 

~ir = ~r{i + z} (2 - 27) 

with u = — s — i. We get: 

a el (S) « a^ A = [ {l ' m '^ S? dx f dhei(x) (2.28) 

Jx m i n Jm%—s dt 

where x = s/S and j e i( x ) is the elastic contribution to the equivalent photon distribution 
of the proton 



- |,J| ~2^ X L» Tl lxix ~ ) + tJ l(t) + 2<t) /' < 2 - 29 ' 



with 



2 2 
m x 



tmin ~ 0° ^-max ~ -, ■ (2.30) 

1 — X 

To clarify the physical meaning of x, let us introduce the variable x 7 : 

x 7 = (2.31) 

It is possible to show that x 7 represents the fraction of the longitudinal momentum of the 
proton carried by the virtual photon, so that one can write 



k = x 7 P + k, (2.32) 



with k ■ P = 0. Using ()2.2|) one gets 



o — m 



which reduces to x in the EPA limit. One can also define the leptonic variable xf. 

»' = Wlh-y (2 - 34) 



where Qf = —i. When t ~ 0, it turns out that also xi 



x. 



III. INELASTIC QED COMPTON SCATTERING 

To calculate the inelastic QED Compton scattering cross section, we extend the approach 
discussed in the previous section. In this case, an electron and a photon are produced in the 
final state with a general hadronic system X. In other words, we consider the process 

e(0 + p(P) - e(Z') + 7 (*') + X{P X ), (3.1) 

where Px = J2xi Px d is the sum over all momenta of the produced hadronic system. Let 
the invariant mass of the produced hadronic state X to be W. Eq. (|2.2|) still holds with 
Q 2 = —t. The cross section for inelastic scattering will be 

<w(S) = — ]_ J dPS 2+N {l + P; I', k', P Xl , Px N )\M incl \ 2 , (3.2) 

where 



I M inel | 2 = ^HZi(P, Px)T^(l, k; I', k') (3.3) 
is the squared matrix element and 

HZi(P,Px) = ~ E E (Px I ^(0) I PT(Px I no) I P) . (3.4) 

spins X 

If we rearrange the (2 + iV)-particle space phase into a sequence of a 2-particle and a iV- 
particle one, we get 

1 r dW 2 r ds r 1 

■US) - ^^) / IT / 27 / + P; l ' + P ^ W ^ k ^ V> 

(3.5) 
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where is given by Eq. (j2.11j) and W^ u is the hadronic tensor for inelastic scattering 



W^ = JdPS N (P-k;P Xl ,....,P XN )H! 



inel 



(3.6) 



The hadronic tensor is parametrized as 



Aire' 



- (QV" + k»k u )F 1 (x B , Q 2 ) + (2x B P» - W)(2x B P v - k l 



F 2 (x B ,Q 2 



2x B 



(3.7) 



where x B is the Bjorken variable given by 

Q 2 



X B 



Q 2 



2P ■ k Q 2 + W 2 -m 2 



(3.8) 



Using 



dPS 2 {l + P;l' + k',P x ) 



dQ 2 



1tt(S — m?) 



as before, we get 



a 



Att(S - m 2 ) 2 Jwli 
S-m 2 



dW z 



(Vs-w) 2 



ds 



> 2 

- max 



d Q 2 ^ /"* 



(3.9) 



S-m 2 
s + Q 2 



x 1 



s + Q 2 



+ {W 2 -m 



, ( i 2 (S - m 2 ) 1 m 2 - W 2 
Q 2 {s + Q 2 )~Q 2 + 2Q 4 



x [3X 1 (s,Q 2 ,t)+X 2 (s,Q 2 ,t)] + ^— (W 2 -m 2 ) + 
x [X^S, Q 2 , t) + X 2 (s, Q 2 , £)] - X 1 (s, Q 2 , t) 
-X 2 (s,Q 2 ,t)F 1 (x B ,Q 2 



(W 2 - m 2 ) 2 2m 



+ 



2Q 4 Q 



F 2 (x B , Q 



2 ^ X B 



(3.10) 

Here Xi(s,Q 2 ,i), with z = 1,2, are given by Eqs. (j2.15j) - (j2.16|) with t replaced by — Q 2 . 
The limits of the integration over Q 2 are: 



Qn„ = -m 2 -W 2 + ^[(S + m 2 )(S-s + W 2 ) ? (S - m 2 ) y/(S - s + W 2 ) 2 - ASW 2 



(3.11) 



while the extrema of t are the same as Eq. (|2.18|) . The limits W^ linmax are given by: 



= (m + m w ) 2 , 



'Vs 



(3.12) 



where m^ is the mass of the pion. F\ and F 2 are the usual structure functions of the proton. 
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In the EPA, we neglect m 2 compared to S and Q 2 compared to s as before. Using Eqs. 
<^2B and (HHBJ) , we can write 

epa /-a-Wv^) 2 ro ^d&fxS,^ /nin , 

cr in ei(<S) ~ cr inel = / diC / rft 7 ine l (x, x6 ) , (3.13) 

Jx min Jm\ — s at 

where again x = s/S and (x, xS) is the inelastic contribution to the equivalent photon 
distribution of the proton [17]: 

a /-i fQ™*dQ 2 yr ^ 1 n -^ 2 ^2^2, 



7kei(s, artf) = [ dy [ 



y J\ y 2 y 2 Q 2 



Following 



(3.14) 

\y / j 

j], [nj] we shall use the LO Callan-Gross relation 

F L {x B ,Q 2 ) = F 2 (x B ,Q 2 )-2x B F 1 (x B ,Q 2 ) = (3.15) 
in our numerical calculations. The limits of the Q 2 integration can be approximated as 



2 2 

2 ^ ^ 2 / \ 

Qmin = 1 5 Qmax = S - (3.16) 
I — X 



Our expression of •yi ne \(x,xS) differs from 11] by a (negligible) term proportional to m 2 . 

IV. NUMERICAL RESULTS 

In this section, we present an estimate of the cross section, calculated both exactly and 
using the equivalent photon approximation of theproton. We have used the same kinematical 
cuts as used by the HI collaboration at HERA [9] , which are slightly different from the ones 
in Q|. They are imposed on the following lab frame variables: energy of the final electron 
E' e , energy of the final photon E' y , polar angles of the outgoing electron and photon, 9 e 
and (9 7 respectively, and acoplanarity angle <f), which is defined as <fi = | 7r — |0 7 — <p e \ |, 
where 7 and <p e are the azimuthal angles of the outgoing photon and electron respectively 
(0 < <py, 4> e < 2 it). The cuts are given by: 

4 GeV, E' e + E! y >20GeV, (4.1) 
0.06 < # e ,# 7 < 7T- 0.06, (4.2) 
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< < \- (4.3) 

The energies of the incoming particles are: E e = 27.5 GeV (electron) and E p = 820 GeV 
(proton). In our conventions, we fix the lab frame such that <p e = 0, so the acoplanarity 
will be <p — | 7r — These cuts reflect experimental acceptance constraints as well as 
the reduction of the background events due to emitted photons with (V + k') 2 ~ and/or 
(1 — k') 2 « 0, which are unrelated to the QED Compton scattering process (for which 
—k 2 = Q 2 « but with both (/' + k') 2 and (/ — k') 2 finite), i.e. photons emitted parallel to 
the ingoing (outgoing) electron or from the hadron vertex . 

We numerically integrate the elastic and inelastic cross sections given by Eqs. ()2.23|) and 
()3.10|) . To implement the cuts in Eqs. ()4.1|) - (j4.3|) . we express E' e , E' , cos9 e , cos # 7 and 
cos0 in terms of our integration variables s, t, i, (p* (and W 2 in the inelastic channel), as 
explained in Appendices A-D. More explicitly, we use Eqs. flB7j) - (jBll|) . ()A7j) . (|A8|) . together 
with Eqs. flAllJ) - (jA13|) for the elastic channel and Eqs. (jC3|) - ()C5|) for the inelastic one. The 
cuts imposed on the lab frame variables restrict the range of our integrations numerically. 
In this way, we are able to remove the contributions from outside the considered kinematical 
region. 

In the calculation of the inelastic cross section, we have used the ALLM97 parametriza- 
tion of E2(xb,Q 2 ) [16], which is obtained by fitting DIS data of HERA and fixed target 
experiments together with the total pp and jp cross sections measured and is expected to 
hold over the entire available range of xb and Q 2 . We have not considered the resonance con- 
tribution separately but, using the so called local duality we have extended the ALLM97 
parametrization of F 2 from the continuous (W > 1.8 GeV) down to the resonance domain 
(m + m-x < W < 1.8 GeV). In this way it is possible to agree with the experimental data 
averaged over each resonance, as pointed out in [9J]. The elastic contribution to the EPA was 
calculated using Eq. ()2.28j) subject to the additional kinematical restrictions given b\ 
(j4.1M.2|) . It corresponds to the same equivalent photon distribution as presented in 
For the inelastic channel we have used Eq. f)3. 13j) together with Eq. ()3.14|) . the cuts being 
the same as in the elastic case. We have taken Fl = and used the ALLM parametrization 
of Fo^_m order to compare consistently with the 'exact' cross section. We point out that in 
Jl F 2 (xb,Q 2 ) in 7i ne ] (x, xS) was expressed in terms of parton distributions for which 
the LO GRV parametrization was used, together with Q min = 0.26 GeV 2 so as to guar- 



M. 
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antee the applicability of perturbative QCD Q]. y me i(x,xS) in our case gives slightly higher 
results than the ones obtained with the photon distribution presented in B, 0] . 

The Compton process turns out to be dominated by the elastic channel, in fact after 
Monte Carlo integration, we find that cr c i = 1.7346 pb, while <r inel = 1.1719 pb. The 
approximated calculation gives the results: cr^ PA = 1.7296 pb and crP^ = 1.5969 pb. This 
means that in the kinematical region under consideration, the total (elastic + inelastic) 
cross section calculated using the EPA agrees with the exact one within 14% and that the 
approximation turns out to be particularly good in describing the elastic process, for which 
the agreement is within 0.3%. This is not surprising since in the EPA one assumes Q 2 = 0, 
which is not true especially in the inelastic channel and the inelastic cross section receives 
substantial contribution from the non-zero Q 2 region. In terms of the kinematical cuts, 
the EPA corresponds to the situation when the outgoing electron and the final photon are 
observed under large polar angles and almost opposite to each other in azimuth, so that 
the acoplanarity is approximately zero. For elastic scattering there is a sharp peak of the 
exact cross section for = 0, contributions from non-zero <fi are very small in this case. But 
the inelastic cross section receives contribution even from non-zero 0, so that in this case 
the discrepancy from the approximated result is higher. The discrepancy of the total cross 
section with the approximate one is thus entirely due to the inelastic part. 

In Fig. 2 we have compared the total cross sections (exact and EPA) in different x\ bins, 
in the region 1.78 x 10~ 5 < xi < 1.78 x 10 _1 . Fig. 3 shows that the agreement improves 
slightly for bins in the variable Xy. Since Xsy — %i for Q 2 ~ 0, the elastic process is not 
sensitive to this change of variables. We point out again that in the EPA limit (Q 2 = 0) 

In Fig. 4 we show the exact and the EPA cross section in x\ and Qf bins together with the 
experimental results and the estimates of the Compton event generator, already presented 
in jjj. Except for three bins, our exact result agrees with the experiment within the error 
bars. The slight difference of our exact result and the one of |9] may be due to the fact 
that in Si the cross section is calculated using a Monte Carlo generator in a step by step 



iteration 



13j which starts by assuming Q 2 = 0, while we did not use any approximation. 
Our exact result is closer to the EPA in most of the kinematical bins as compared to 
The total cross section in the EPA lies above the 'exact' one in most of the bins. 

For completeness, we have shown the numerical values of the exact and EPA double 
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differential cross sections, both for the elastic (Table 1) and inelastic (Table 2) contributions. 



The kinematical bins are the same as in [9( . The exact results when the bins are in x 7 instead 
of Xi are also shown. The EPA elastic cross section agrees within 1% with the exact one 
for all the x\ bins. The agreement becomes slightly better if we consider x 7 bins. For 
the inelastic channel, the discrepancies from the EPA results are considerably higher. Our 
'exact' results lie closer to the EPA compared to in almost all the bins. The result in x 7 
bins show better agreement with the EPA compared to the x\ bins, especially for higher x 7 . 
The discrepancy with the EPA is about 20 — 30% in most of the bins, higher in some cases. 



V. SUMMARY 



In this work, we have calculated both elastic and inelastic QED Compton scattering cross 
section in the unpolarized case. Our approach for the total cross section is manifestly co- 
variant and we have used the same cuts as in the HERA-HI experiment. The numerical 
estimates of the exact cross section for different kinematical bins are presented and compared 
with the EPA and the experimental results. The exact cross section in the elastic channel 
agrees within 1% with the approximate one. The discrepancy is thus due to the inelastic 
channel. The discrepancy with the Monte Carlo estimate of P| is also shown. For both 
elastic and inelastic cross sections, our exact result is closer to the EPA as compared to |^|. 
The agreement is even better if the bins are in x 7 instead of x\. Our approach can be ex- 
tended to calculate the corresponding cross section for polarized scattering and also for other 
polarized and unpolarized processes having photon induced subprocesses, in order to check 
how accurately the cross section is given by the equivalent photon approximation. Also this 
would predict the kinematical cuts necessary for the extraction of 7(2;, Q 2 ) experimentally 
from these processes. 
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APPENDIX A: KINEMATICS IN THE e - 7 C. M. FRAME (ELASTIC) 

In this appendix, we discuss the kinematics of the QED Compton scattering in the c. m. 
frame of the outgoing e — 7 system. The 4-momenta of the particles are given by: 
incident electron: I = (E*, 0, 0, — E*), 

incident proton: P = (E*, P*sin9*, 0, P*cos9*), where P* = \J E* 2 - m 2 , 

outgoing electron: I' = E'*(l, sin 9* cos tp* , sin 9* sin tp* , cos 6*), 

outgoing photon: k' = E'*(l, — sin 0* cos tp* , — sin 9* sin tp* , — cos 8*). 

We have also: 4-momentum of the virtual photon: k = (El, 0, 0, E* e ) , with k 2 = t and 

E* k = y/Ef + t. (Al) 

The overall momentum conservation allows us to write the 4-momentum of the final proton 

as 

P' = l + P-l' -k'. (A2) 
We introduce the following Lorentz invariants 

s = (l' + k') 2 = AE'* 2 , (A3) 

t = (l- I') 2 = -2E* e E'*(l + cose*), (A4) 

u — (I — k') 2 = -2E*E'*{l-cos9*), (A5) 

S = (I + Pf = m 2 + 2E* p E* e + 2E* e P; cos 9* p , (A6) 

T = (P -I') 2 = m 2 - 2E'*(E; - P* sin 9* sin 9* p cos tp* - P p * cos 9* cos 9*), (A7) 

U =(P- k') 2 = m 2 - 2E'*(E; + P* sin 9* sin 9* p cos tp* + P p * cos 9* cos 9* p ). (A8) 
In addition they satisfy: 

§ + t + u = t, S + T + U = -t + 3m 2 . (A9) 
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Using the relations above, it is possible to write the energies of the particles in the lab frame 
in terms of the integration variables s, t , t and the constant S: 

K = El = (A10) 

and 



„ S-m> + t „, y/(S - m? + if - 4sm* ... 
^ = ^|-' P;= Wl ' E (AU) 



Similarly, for the angles we have 



cosO' = f s "_ t 2t , (A12) 



P (S - t)[(S + t - m 2 ) 2 - 4sm 2 ]^ 

In particular Eqs. (jA7|) and ()A8|) . through Eqs. ()A11|) - ()A13|) . express T and U in terms of 
our integrations variables s, t, t, (p* and we have used them to relate the lab frame variables 
to the integration ones, as shown in the next appendix. 

APPENDIX B: KINEMATICS IN THE LAB FRAME (ELASTIC) 

In this appendix, we give the kinematics of the same process in the lab frame. The 4 
momenta of the particles are given by: 
incident electron: I = (E e , 0, 0, —E e ), 



incident proton: P = (E p , 0, 0, P p ), where P* = y E p 2 — m 2 , 

outgoing electron: I' = E' e (l, sin^e, 0, cos# e ), 

outgoing photon: k! = E' (1, sin# 7 cos0 7 , sin# 7 sin0 7 , cos# 7 ). 

Here we have chosen the frame such that the outgoing electron has zero azimuthal angle. 
The Lorentz invariants are: 

s=(l' + k'f = 2E' e E! y (l - sinfl e sin# 7 cos0 7 - cos# e cos# 7 ), (Bl) 
t = (l-l') 2 = -2E e E' e (l + cos6 e ) } (B2) 



u — {I — k') = —1E e E'{\ + cos# 7 ), (B3) 
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S = (I + Pf = m 2 + 2E e (E p + P p ), 



(B4) 



T = (P- I') 2 = m 2 - 2E' e (E p - P p cos6 e ), (B5) 



U = {P- k') = m 2 - 2E' 1 {E P - P p cosfl 7 ). (B6) 

The polar angles in the lab frame can be written in terms of the invariants and the incident 
energies: 

E p t-E e (T-m 2 ) 

co ^ = p ri+E .(T- m >y (B7) 

7 P p (t-s-t) + E e (U - m 2 ) v ' 

In the same way, for the energies of the final electron and photon we have: 

= _ iPr + E , {T 

O — 771 

1 S — m 2 

The azimuthal angle of the outgoing photon is: 

2^;(l-cosg e cosg 7 )-s 

cos©-, = ■ , (Bll) 

Vl 2E' e E' 1 sm6 e sm6 1 ' v ; 

which is related to the acoplanarity angle by the relation = | n — 7 | . Using Eqs. (|A7jl 
and JSHj) together with Eqs. (|A11J) - (|A13J) . the formulae above for cos^e, cos6* 7 , -Eg , E'^ and 
cos0 7 can be expressed in terms of the invariants and the incident energies. Eqs. (|B7|) - (jBll|) 
are needed to implement numerically the kinematical region under study, since they relate 
the lab variables (energies and angles) to the ones used for the integration. 



APPENDIX C: KINEMATICS IN THE e — 7 C. M. FRAME (INELASTIC) 

In this frame most of the expressions remain the same as given in appendix A. The 
relations among the invariants are now: 

O 2 

§ + t + u = -Q 2 , S + T + U = 3m 2 + ^L, (CI) 

x B 
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where Q 2 = —t and xb can be written 



as 



XB = TV) i wo 2" ( C2 ) 

Q 2 + W 72 — m 2 

The only formulae which are different are the ones involving E* and cos#*. So Eq. (jAll|) 
will be replaced by 



p* g-Q 2 -^ 2 p * y/(g-g 2 -^ 2 ) 2 -4gm 2 VI rp ,. 

Ep = Wl ' P= Ws ' (C3) 



while for the angles: 



q 2 + s + 2t 

COs9 =- s + Q 2 ' (C4) 



and 

fl * 2g(g-m 2 )-(g-g 2 -^ 2 )(g + Q 2 ) 
cosfcL = j — . (Co) 

(s + Q 2 )[(S - Q 2 - W 2 ) 2 - Asm 2 }* 

Eqs. O-fJCSJ) reduce to the Eqs. (fATTl)-(fA"i"5|) of the elastic channel for W = m and 
Q 2 = -t. 



APPENDIX D: KINEMATICS IN THE LAB FRAME (INELASTIC) 

The invariants in the case of inelastic scattering are the same as in the elastic case: Eqs. 
(|Blj) - ([B6p . Eq. (jCij) describes the relation among them. The expressions of cos^g, cos# 7 , 
E' e , E' and cos</> 7 , in terms of the integration variables W 2 , s, Q 2 , t, <p* are given by Eqs. 
(Ib7|) - (IbTT1) together with Eqs. (|A"7|l . (|A"8j) as before, but now Eqs. (fH3 |) -(|H5 |) will replace 
Eqs. (lATH - (lAT3l) . 
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Fig. 2: Cross section for Compton process at HERA-HI. The cuts applied are as described 
in the text. 
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Fig. 3: Cross section for Compton process at HERA-HI. The bins are in x 1 . The cuts 
applied are as described in the text. 
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Fig. 4: Double differential cross section for QED Compton scattering at HERA-HI. The 
data are taken from The kinematical bins correspond to Table 1. The continuous line 
corresponds to our exact calculation, the dotted line to the calculation in yj and the dashed 
line to the EPA. 
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x bin 


Qf bin 






a e\ 


..EPA 
el 


1.78 x 1(T 5 - 5.62 x 1(T 5 
1.78 x 1CT 5 - 5.62 x 1(T 5 


1.5 - 2.5 
2.5 - 3.5 


2.428 x 10 2 
5.099 x 10 1 


2.342 x 10 2 
4.71 x 10 1 


2.446 x 10 2 
5.201 x 10 1 


2.461 x 10 2 
5.051 x 10 1 


5.62 x 1(T 5 - 1.78 x 1(T 4 
5.62 x 1(T 5 - 1.78 x 10~ 4 
5.62 x 10" 5 - 1.78 x 10" 4 


1.5 - 5.0 
5.0-8.5 
8.5 - 12.0 


5.279 x 10 2 
2.396 x 10 2 
8.496 x 10 1 


5.319 x 10 2 
2.327 x 10 2 
8.32 x 10 1 


5.259 x 10 2 
2.404 x 10 2 
8.559 x 10 1 


5.247 x 10 2 
2.395 x 10 2 
8.571 x 10 1 


1.78 x 10" 4 - 5.62 x 1(T 4 
1.78 x 1(T 4 - 5.62 x 10~ 4 
1.78 x 10" 4 - 5.62 x 10" 4 


3.0 - 14.67 
14.67 - 26.33 
26.33 - 38.0 


2.080 x 10 2 
1.373 x 10 2 
3.712 x 10 1 


2.036 x 10 2 
1.388 x 10 2 
3.86 x 10 1 


2.056 x 10 2 
1.373 x 10 2 
3.720 x 10 1 


2.061 x 10 2 
1.372 x 10 2 
3.695 x 10 1 


5.62 x 10" 4 - 1.78 x 10~ 3 
5.62 x 1(T 4 - 1.78 x 1(T 3 
5.62 x 10~ 4 - 1.78 x 10~ 3 


10.0 - 48.33 
48.33 - 86.67 
86.67 - 125.0 


5.947 x 10 1 
3.714 x 10 1 
1.056 x 10 1 


5.71 x 10 1 
3.85 x 10 1 
1.028 x 10 1 


5.918 x 10 1 
3.715 x 10 1 
1.057 x 10 1 


5.921 x 10 1 
3.704 x 10 1 
1.054 x 10 1 


1.78 x 1(T 3 - 5.62 x 1(T 3 
1.78 x 10~ 3 - 5.62 x 1(T 3 
1.78 x 1CT 3 - 5.62 x 10 -3 


22 - 168 
168 - 314 
314 - 460 


1.913 x 10 1 
1.239 x 10 1 
5.917 


1.877 x 10 1 
1.229 x 10 1 
6.02 


1.909 x 10 1 
1.239 x 10 1 
5.915 


1.909 x 10 1 
1.238 x 10 1 
5.914 


5.62 x 1(T 3 - 1.78 x 1(T 2 
5.62 x 10~ 3 - 1.78 x 1(T 2 
5.62 x 1CT 3 - 1.78 x 1(T 2 


- 500 
500 - 1000 
1000 - 1500 


4.811 
9.271 
2.572 


5.76 
9.22 
2.65 


4.890 
9.264 
2.571 


4.890 
9.271 
2.573 


1.78 x 10~ 2 - 5.62 x 1(T 2 
1.78 x 1(T 2 - 5.62 x 10~ 2 
1.78 x 1(T 2 - 5.62 x 1(T 2 


- 1500 
1500 - 3000 
3000 - 4500 


8.238 x 10" 1 

2.431 
6.336 x lO" 1 


6.8 x 10" 1 

2.69 
7.7 x 10- 1 


9.085 x lO" 1 

2.430 
6.328 x 10" 1 


9.086 x lO" 1 

2.434 
6.345 x 10" 1 


5.62 x 10~ 2 - 1.78 x 1CT 1 
5.62 x 1(T 2 - 1.78 x KT 1 
5.62 x 10~ 2 - 1.78 x 1CT 1 


10 - 6005 
6005 - 12000 
12000 - 17995 


3.120 x 10 _1 
2.437 x 10" 1 
0.000 


4.27 x 10 _1 
2.13 x 10" 1 
0.000 


3.120 x 10 _1 
2.438 x 10" 1 
0.000 


3.117 x 10 _1 
2.436 x 10" 1 
2.461 x 10~ 2 



TABLE I: Double differential (elastic) QED Compton scattering cross section. a e \ is the exact 
result in Eq. ()2.23|) . (7^f n corresponds to the results in 9]. The x-bins refer to x\ in Eq. (|2.34l) 
except for cr*j where they refer to x 7 in Eq. (|2.33f> . c^ PA is given in Eq. (|2.28f) where 
is expressed in GeV 2 and the cross-sections are in pb. 
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x bin 


Qf bin 


Cinel 




* 

inel 


O ; ne [ 


1.78 x 1(T 5 - 5.62 x 1(T 5 
1.78 x 1CT 5 - 5.62 x 1(T 5 


1.5 - 2.5 
2.5 - 3.5 


7.996 X 10 1 
2.880 x 10 1 


6.170 x 10 1 
2.050 x 10 1 


7.503 x 10 1 
4.142 x 10 1 


1.529 x 10 2 
3.361 x 10 1 


5.62 x 1(T 5 - 1.78 x 1(T 4 
5.62 x 1(T 5 - 1.78 x 10~ 4 
5.62 x 10" 5 - 1.78 x 10" 4 


1.5 - 5.0 
5.0-8.5 
8.5 - 12.0 


2.361 x 10 2 
1.139 x 10 2 
3.442 x 10 1 


2.296 x 10 2 
1.062 x 10 2 
2.890 x 10 1 


2.364 x 10 2 
1.500 x 10 2 
6.291 x 10 1 


4.116 x 10 2 
2.099 x 10 2 
8.094 x 10 1 


1.78 x 10" 4 - 5.62 x 1(T 4 
1.78 x 1(T 4 - 5.62 x 10~ 4 
1.78 x 10" 4 - 5.62 x 10" 4 


3.0 - 14.67 
14.67 - 26.33 
26.33 - 38.0 


2.119 x 10 2 
1.005 x 10 2 
2.868 x 10 1 


2.287 x 10 2 
9.230 x 10 1 
2.570 x 10 1 


1.278 x 10 2 
1.220 x 10 2 
4.682 x 10 1 


2.210 x 10 2 
1.556 x 10 2 
4.558 x 10 1 


5.62 x 10" 4 - 1.78 x 10~ 3 
5.62 x 1(T 4 - 1.78 x 1(T 3 
5.62 x 10~ 4 - 1.78 x 10~ 3 


10.0 - 48.33 
48.33 - 86.67 
86.67 - 125.0 


1.079 x 10 2 
4.779 x 10 1 
1.315 x 10 1 


1.064 x 10 2 
4.590 x 10 1 
1.132 x 10 1 


5.383 x 10 1 
5.690 x 10 1 
2.050 x 10 1 


8.092 x 10 1 
5.272 x 10 1 
1.587 x 10 1 


1.78 x 1(T 3 - 5.62 x 1(T 3 
1.78 x 10~ 3 - 5.62 x 1(T 3 
1.78 x 1CT 3 - 5.62 x 10 -3 


22 - 168 
168 - 314 
314 - 460 


4.758 x 10 1 
2.010 x 10 1 
6.940 


4.917 x 10 1 
1.735 x 10 1 
5.760 


2.426 x 10 1 
2.350 x 10 1 
9.303 


3.080 x 10 1 
2.058 x 10 1 
1.021 x 10 1 


5.62 x 1(T 3 - 1.78 x 1(T 2 
5.62 x 10~ 3 - 1.78 x 1(T 2 
5.62 x 1CT 3 - 1.78 x 1(T 2 


- 500 
500 - 1000 
1000 - 1500 


1.482 x 10 1 
1.228 x 10 1 
3.135 


1.432 x 10 1 
9.890 
2.600 


7.327 
1.318 x 10 1 
3.756 


8.825 
1.687 x 10 1 
4.885 


1.78 x 10~ 2 - 5.62 x 1(T 2 
1.78 x 1(T 2 - 5.62 x 10~ 2 
1.78 x 1(T 2 - 5.62 x 1(T 2 


- 1500 
1500 - 3000 
3000 - 4500 


3.585 
3.778 
9.675 x lO" 1 


2.500 
2.150 
6.600 x 10" 1 


1.623 
3.866 
1.088 


1.811 
4.867 
1.341 


5.62 x 10~ 2 - 1.78 x 1CT 1 
5.62 x 1(T 2 - 1.78 x KT 1 
5.62 x 10~ 2 - 1.78 x 1CT 1 


10 - 6005 
6005 - 12000 
12000 - 17995 


1.093 
5.496 x 10" 1 
6.164 x 10~ 2 


1.460 x 10 _1 
2.110 x 10" 1 
4.300 x 10~ 2 


6.589 x 10 _1 
5.782 x lO" 1 
0.000 


7.147 x 10 _1 
5.922 x 10" 1 
6.791 x 10~ 2 



TABLE II: Double differential (inelastic) QED Compton scattering cross section. <7i ne i is the exact 
result in Eq. (|3.10|) . corresponds to the results in The x-bins are as in Table I, i.e. refer 
to x\ in Eq. Q2.34JI except for af nel where they refer to x 7 in Eq. (j2.33f) . crP^j is given in Eq. (|3.13|) 
where x = x~. Qf is expressed in GeV 2 and the cross-sections are in pb. 



